Study of phase diagram and superconducting states in LaFeAsOi- 
based on the multiorbital extended Hubbard model 



hHt 



cn 

o 

<N 

< 



o 
o 

u 

Oh- 



o 



> 

m 
m 
0\ 
cn 

o 
m 



^ 



Y. Yamakawa^, S. Onari^, H. Kontani^, N. Fujiwara^, S. Iimura^, and H. HosoNO^ 

^ Department of Physics, Nagoya University, Furo-cho, Nagoya 464-8602, Japan. 

^ Department of Applied Physics, Nagoya University, Furo-cho, Nagoya 464-8602, Japan. 

^ Graduate School of Human and Environmental Studies, Kyoto University, 

Yoshida-Nihonmatsu-cho, Sakyo-ku, Kyoto 606-8501, Japan. 

^ Materials and Structures Laboratory, Tokyo Institute of Technology, 

4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan. 

(Dated: April 16, 2013) 

To understand the recently established unique magnetic and superconducting phase diagram of 
LaFeAsOi-ccHa;, we analyze the realistic multiorbital tight-binding model for x = ^ 0.4 beyond the 
rigid band approximation. Both the spin and orbital susceptibilities are calculated in the presence of 
the Coulomb and charge quadrupole interactions. It is found that both orbital and spin fluctuations 
strongly develop at both x ^ ^ and 0.4, due to the strong violation of the rigid band picture 
in LaFeAsOi-ccHcc. Based on this result, we discuss the experimental phase diagram, especially 
the double-dome superconducting phase. Moreover, we show that the quadrupole interaction is 
effectively produced by the vertex correction due to Coulomb interaction, resulting in the mutual 
development of spin and orbital fluctuations. 

PACS numbers: 74.70. Xa, 74.20.-z, 74.25. Dw 



Since the discovery of high-Tc superconductivity in 
Fe-based superconductors [1], its pairing mechanism has 
been studied very intensively. In unconventional super- 
conductors, the phase-diagram in the normal state gives 
us very important hints to understand the mechanism 
of superconductivity. In many heavy fermion supercon- 
ductors, for example, the superconducting (SC) phase is 
next to the magnetic ordered phase, indicating the occur- 
rence of spin-fluctuation mediated superconductivity. In 
Fe-based superconductors, in contrast, the ferro-orbital 
order occurs in the orthorhombic phase ^2|, [3|, and the 
structure or orbital instabilities are realized in the nor- 
mal state 0, [5| , in addition to the magnetic instability. 
Based on this fact, both the spin- fluctuation mediated 
5±-wave state [6-9] and the orbital-fluctuation mediated 
5++-wave state [1Q|-[12] have been discussed. The former 
(latter) SC gap with (without) sign reversal is induced 
by repulsive (attractive) interaction between electron-like 
and hole- like Fermi surfaces (FSs). 

Recently, very unique phase diagram of H-doped 
LaFeAsO, LaFeAsOi-a^H^^, is determined in Ref. [13|: 
The structure and magnetic transitions in mother com- 
pound are replaced with the SC phase at x '^ 0.03, and 
interesting double-dome structure of T^ is obtained be- 
tween X = 0.03 and < 0.5. The maximum Tc of the 
first (second) dome is about 25 K (40 K) at x ~ 0.1 
{x r^ 0.35). For x > 0.4, recent NMR measurement [ij 
detected the incommensurate magnetic order, in addition 
to the highly anisotropic electric field gradient that indi- 
cates the occurrence of the non-magnetic orbital order. 

In LaFeAsO i_aj Ha;, the electron filling per Fe is n = 
6-\-x since each H-dopant becomes H~ ion, and therefore 
the electronic states are expected to be similar to those of 
LaFeAsOi-a^Fa; [l3[. Based on this fact, the band struc- 
ture and FSs for x = ^ 0.4 had been derived from the 



local-density-approximation (LDA) band calculation us- 
ing the virtual crystal approximation in Ref. [13|. It is 
found that the rigid band picture is no more vaHd, since 
the band structure is strongly modified with increasing 
F-doping. The derived realistic band structure now en- 
ables us to perform quantitative theoretical study of the 
pairing mechanism of 1111 systems. 

In this paper, we study the electronic and SC states 
in LaFeAsO 1- a; Ha; for x = ~ 0.4, by constructing 
the realistic multiorbital models beyond the rigid band 
approximation. Using the random-phase-approximation 
(RPA), both the spin and orbital susceptibilities are cal- 
culated in the presence of the Coulomb interaction U 
and charge quadrupole interaction g. Assuming mono- 
tonic x-dependencies of these interactions, strong spin 
and orbital fluctuations are obtained for both x ^ and 
X ~ 0.4. Based on this result, the origin of the double- 
dome structure of Tc in LaFeAsO i-a^H^^ is discussed, by 
applying both the orbital-fluctuation mediated s++-wave 
scenario and spin-fluctuation mediated 5±-wave one. We 
discuss that g is effectively induced by the vertex correc- 
tion (VC) of the Coulomb interaction beyond the RPA. 

In LaFeAsOi-a^Fa;, the spin fluctuations observed by 
NMR and neutron inelastic scattering are very small 
in slightly over-doped compounds {x ^ 0.08). In the 
over-doped compound with x = 0.14, Tc increases from 
20K to 43K by applying 3.7GPa, irrespective that 1/^ 
remains very small independently of the pressure [l5l |. 
These facts indicate the weak correlation between Tc 
and spin fluctuations. Moreover, impurity effect on Tn 
is very small in both 1111 [jj, and 122 [H, [Tf 
compounds, indicating the realization of the s++-wave 
state [20, UM' Unfortunately, these experiments on 
LaFeAsO i_a^ Ha; have not been performed yet. 

First, we perform the LDA band calculation for 
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FIG. 1: (Color online) (a) Band structures for x = (solid 
line) and x — 0.4 (dashed line) in the present model, (b) FSs 
for X = and (c) FSs for x — 0.4. The weights of xz-^ yz- and 
xy-orbital on the Fermi surfaces are represented by diameter 
of green-, red-, and blue-circles, respectively, (d) Weight of 
orbitals on the Fermi surface around X point for x = 0.4 as 
a function of azimuthal angle 0. (e) FSs for n = 6.4 given by 
the rigid band model for LaFeAsO. 



the three xz-^ yz- and x?/-orbitals: The orbital character 
of the electron-pocket for x = 0.4 is shown in Fig. [1] (d). 
In the case of x = 0, the electron-hole (e-h) FS nesting 
with the nesting vector Q = (tt, 0), (0,7r) is the most 
important. On the other hand, in the case of x = 0.4, 
electron-electron (e-e) FS nesting Q ~ (7r,7r/3), (7r/3,7r) 
is more important since the hole pockets become very 
small. In both x = and x = 0.4, both the intra-orbital 
nesting and inter-orbital nesting (between xz/yz and 
xy) are important. Then, the former (latter) nesting 
gives rise to the strong spin (orbital) fluctuations, as 
discussed in Ref. [iQj . 

Fig. [1] (e) shows the FSs for n = 6.4 given by the 
model parameters for x = (LaFeAsO). In this "rigid- 
band approximation", the hole pockets disappear, and 
the e-e FS nesting is worse since the shape of the electron- 
pockets are more rounded. 

Next, we explain the interaction term. We intro- 
duce both the Coulomb interaction (/7, U\ J, J') and 
quadrupole interaction (g). The latter interaction is 

site 

^quad = -g{coi) J2 (oL-oL + oU-oU + oiy-oiy) , 

i 

(1) 

where g{uji) = g-uj'^/{ujf-\-uj'^): g = g{0) is the quadrupole 
interaction at ui = 2l7rT = 0, and Uc is the cutoff en- 
ergy. Or is the quadrupole operator [lOi] , which has many 
non-zero off-diagonal elements. By introducing small 
g {^ 0.2eV), strong Oxz,yz-'^yp^ antiferro-quadrupole 
fluctuations are caused by the good inter-orbital nesting, 
as explained in Ref. [1Q| |. 

Now, we perform the RPA for the present model at 
T = 0.02eV, by assuming that J = J' and [/ = [/' + 2 J, 
and fix the ratio J/U = 1/6. We use 64 x 64/c meshes 
and 512 Matsubara frequencies. We set the unit of energy 
as eV hereafter. The spin (orbital) susceptibility in the 
RPA is given by 



LaFeAsOi_^F^ for x = - 0.4 using WIEN2K code with 
virtual crystal approximation, where the oxygen sites are 
substituted for virtual atoms with a fractional nuclear 
charge. Since H~^ ion acts as F~^ ion, the electronic 
states of LaFeAsO 1- a? Ha; are expected to be similar to 
those of LaFeAsO i_aj Fa; [13]. Next, we derive the two- 
dimensional five-orbital tight-binding models for each x 
using WANNIER90 code and WIEN2WANNIER inter- 
face [22]: ^° = E, 



r^'HQ) = x'{Q)\i-r'^'HoJi)x'{q) 



(2) 



hk'^^ijcj^kma, where /,m = 1 - 5 



klmcT k kla 

represent the d orbitals with the order 3z'^ — r'^^ xz^ yz^ 
xy, and x^ — y^'- Here, we set x and y axes parallel to 
the nearest Fe-Fe bonds. 

Figure [1] (a) shows the band structures of 
LaFeAsO 1- a; Fa; in the present model. It is obvious 
that band structure for x = 0.4 cannot be reproduced 
by the rigid band shift from that for x = 0. The 
corresponding FSs are shown in Fig. [1] (b) and (c). 
Here, Pi and ^2 are the electron-pockets, and ai, 0^2 and 
7 are the hole-pockets, both of which are composed of 



where q = (q^uJi), and f^^^^ is the interaction matrix for 
the spin (charge) channel composed of /7, U', J, J' and 
9{uJi) M- X° (5) = -jj Efe Gi^ik + q) Gru'v (k) is the 
irreducible susceptibility, where G{k) = [icn + /i — hk]~^ 
is the bare Green function, and e^ = (2n + l)7rT. 

The magnetic (orbital) order is realized when the spin 
(charge) Stoner factor Qfs(c) 7 which is the maximum eigen- 
value of f"(^)x^°H<?.0). is unity. In the RPA, the en- 
hancement of x^ is mainly caused by the intra-orbital 
Coulomb interaction /7, using the "intra-orbital nesting" 
of the FSs. On the other hand, the enhancement of x^ 
in the present model is caused by the quadrupole inter- 
action in Eq. ([1]), utilizing the "inter-orbital nesting" of 
the FSs. 

Now, we study the development of the spin suscep- 
tibility by [/, by putting ^ = 0. Figure [2] shows the 
total spin susceptibility x {q^^ = 0) = T.i,mXti^mmi^) 



for X = 0,0.14,0.24,0.4, by choosing the U to reahze 
as = 0.98. For (a) x = 0, X^{q) has commensurate 
peaks at gr = (tt, 0) and (0, tt) due to the e-h FS nesting. 
These peaks change to incommensurate for (b) x = 0.14, 
reflecting the size imbalance between electron- and hole- 
pockets. As increasing the doping further, the e-e FS 
nesting and e-h FS nesting become comparable. Because 
of the fact, x^{q) fo^" (c) x = 0.24 shows the double-peak 
structure. For (d) x = 0.4, x^{q) shows the incommen- 
surate peak structure due to the e-e FS nesting only. 




FIG. 2: (Color online) q-dependence of x^(g) with as — 0.98 
for (a) X = 0, (b) X = 0.14, (c) x = 0.24 and (d) x = 0.4. 

Figure [3] (a) shows the x-dependence of Uc^ which is 
the critical value of U for the spin order given by the con- 
dition 0^5 = 1. We stress that Uc in the present model is 
much smaller than that in the rigid band model, reflect- 
ing the good e-h (e-e) FS nesting for x < 0.24 {x > 0.24) 
in the present model. Moreover, Uc takes the maximum 
value at x ~ 0.1, and monotonically decreases by depart- 
ing from X ~ 0.1. For this reason, we can explain the 
magnetic orders at x ~ and 0.4, by assuming a simple 
monotonic x-dependence of the interaction: Here, we in- 
troduce U{x) by the liner interpolation between Uc for 
X = and that for x = 0.4, as shown in Fig. [3] (a). The 
x-dependence of U{x) might be explained by the change 
in the d-electron Wannier functions [23], or the under- 
estimation of the bandwidth given by the first-principle 
study for x ~ 0.4. 

We also analyze the SC state using the linearized 
Eliashberg equation: 



XE^wik) 



N 



Y.^^rn.^^^l'^^-^')^'m.rn.^k') 



k' ,mi 



X A^,.^„{k')G'(-k') + 5Sf,,(e„),(3) 



where Aiif{k) is the gap function, and Ae is the eigen- 
value that reaches unity at T = Tc. When T is fixed, 
the larger value of Ae would correspond to the higher 
Tc. Here, W is the pairing interaction given by the RPA, 
{G')~^ = {G)~^ — ^S^ is the normal Green function with 
impurity-induced normal self-energy ^S"^, and ^S^ is 
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FIG. 3: (Color online) (a) Uc in present models (solid line) 
and that in the rigid band model (dotted line) against x for 
^ = 0. U{x) is determined by liner interpolation between Uc 
for X = and that for x = 0.4. (b) Obtained A^; for the 
s±-wave state with nimp = 0%, 2%, 4%. The x-dependence of 
as is also shown. 
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FIG. 4: (Color online) (a) Qc in present models (solid line) 
and that in the rigid band model (dotted line) against x for 
[7 = 0. g{x) is determined by liner interpolation between Qc 
for X = and that for x = 0.4. (b) Obtained A^; for the s++- 
wave state with cjc = 0.02 ~ 1 (nimp = 0). The x-dependence 
of ac is also shown. 



the impurity-induced anomalous self-energy. The expres- 
sions are given in Ref. [10|. Hereafter, we put the orbital- 
diagonal on-site impurity potential as / = 1. Figure [3] (b) 
shows the obtained Ae as a function of x. We see that Ae 
has two peaks at both x = 0.04 and x = 0.36 since x^iQ) 
develops divergently toward x ^ and 0.4. Then, the 
double-dome behavior of Tc in LaFeAsOi-a^Ha^ would be 
explained, since Tc at both boundaries will be suppressed 
by the magnetic and orbital orders. Here, the hole- 
pockets still exist even at x = 0.4, and then the e-e FS 
nesting Q ^ (7r,7r/3) contributes to the 5±-wave state. 
The d-wave state is expected if the hole-pockets disap- 



pear [2J-|26|, while the relation AE(5±-wave) > XE{d- 
wave) is realized for x < 0.4 in the present study. In 
both states, Ae is quickly suppressed when the impurity 
concentration nimp is finite, meaning that the s±- and 
d-wdive states are fragile against impurities. 

In the next stage, we study the development of the 
orbital susceptibility by ^, by putting U = 0. Figure [H 
(a) shows the x-dependence of ^c, given by the condition 
Qfc = 1. Similarly to Fig [3] (a), we introduce g{x) by 
liner interpolation between gc for x = and that for 
X = 0.4 in Fig. [H (a). Then, the orbital ordered states 
are realized for both x ^ and x ~ 0.4. The obtained 
X2424(9) t)y the RPA for x = 0.4 is shown in Fig. [5] (a). 
We also study the orbital-fluctuation-mediated 5++-wave 
state for g = g[x): Figure [41(b) shows the obtained Ae as 
a function of x. The behavior of Ae in the 5++-wave state 
is similar to that of the 5±-wave state because orbital 
susceptibilities develop divergently toward x ^ and 
0.4. For this reason, the double-dome behavior of Tc can 
be also explained by the orbital fluctuations. The value 
of Ae increases for larger cutoff energy ujc'. In Ref. flOj . 
we put ujc = 0.02 since we considered the quadrupole 
interaction due to the Fe-ion oscillations. However, ooc 
for the effective quadrupole interaction due to VC |12| 
depends on the electronic state and not unique. In both 
cases, we should use larger ujc since the used temperature 
(T = 0.02) is much higher than the real Tc. 



(a) /2424(^) X=0.4 
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FIG. 5: (Color online) (a) X2424(9) given by the RPA us- 
ing the quadrupole interaction in Eq. ([T]) for ac = 0.98 (b) 
X2A2a{q) given by the SC-VC method using the Hubbard in- 
teraction for ac — Ois — 0.98. 

Finally, we explain that the spin and orbital fluctua- 
tions mutually develop by taking the VC into account, 
as actually observed in various Fe-based superconduc- 
tors 0, Q. Beyond the RPA, x^{q) in Eq. ^ is re- 



placed with x°(g) + X^(^)(g), where X^^^) is the VC due 
to the Coulomb interaction for the spin (charge) chan- 
nel. In the RPA, the VC is neglected irrespective of its 
importance. We developed the self-consistent VC (SC- 
VC) method in Ref .[12], and found that the Aslamazov- 
Larkin-type VC gives the large mode-coupling between 
spin and orbital, and the spin fluctuations induce strong 
orbital fluctuations due to this coupling. Figure [5] (b) 
shows the strong development of antiferro-orbital fluctu- 
ations, X2424(9) ^ I5 given by the SC-VC method for 
X = 0.4, [7 = 1.1, J/U = 0.073 and ^ = 0. Therefore, 
the RPA analysis using U{x) and g{x) is justifled by the 
SC-VC theory. 

In addition, strong ferro-orbital fluctuations, that cor- 
respond to the orthorhombic structure transition, are 
also induced by the Aslamazov-Larkin-type VC for x ~ 
[12I, [23 . The orbital fluctuations due to VC also occur 
in a simple two-orbital model, which would explain the 
nematic order in Sr3Ru207 [28, 29]. We have recently 
developed the SC-VCH method, in which both the self- 
energy and the VC are taken into account [30]. Using 
this method, the 5++-wave state can be realized for re- 
alistic parameters {J/U ~ 0.1) for x ^ even for ^ = 0. 
It is an important future problem to analyze the present 
model using the SC-VCS method. 

In summary, we have explained the reappearance of 
the spin and orbital orders in LaFeAsOi-a^Ha^ at x ~ 
and X ~ 0.4. Both spin and orbital orders originate from 
the commensurate e-h FS nesting for x ~ 0, and incom- 
mensurate e-e FS nesting for x ~ 0.4. Due to strong 
spin and orbital fluctuations at x ~ and 0.4, both 
the spin-fluctuation mediated 5±-wave state and orbital- 
fluctuation mediated 5++-wave state can be realized, de- 
pending on the magnitude relation of these fluctuations. 
Since small impurity effect on Tc for the flrst SC dome 
\m . flTj indicates the 5++-wave state, the 5++ -^ s± 
crossover will occur with dopin g, in case that the second 
SC dome is the s±-wave state [10]. Thus, the impurity 
effect study for x ~ 0.4 is highly required. 

This study has been supported by Grants-in-Aid for 
Scientiflc Research from MEXT of Japan. The part of 
Tokyo Tech was supported by JSPS First Program. Nu- 
merical calculations were partially performed using the 
Yukawa Institute Computer Facility. 



[1] Y. Kama, T. Watanabe, M. Hirano, and H. Hosono, J. 

Am. Chem. Soc. 130, 3296 (2008). 
[2] T. Shimojima, K. Ishizaka, Y. Ishida, N. Katayama, K. 

Ohgushi, T. Kiss, M. Okawa, T. Togashi, X.-Y. Wang, 

C.-T. Chen, S. Watanabe, R. Kadota, T. Oguchi, A. 

Chainani, and S. Shin, Phys. Rev. Lett. 104, 057002 

(2010). 
[3] M. Yi, D. H. Lu, J.-H. Chu, J. G. Analytis, A. P. Sorini, 

A. F. Kemper, S.-K. Mo, R. G. Moore, M. Hashimoto, 

W. S. Lee, Z. Hussain, T. P. Devereaux, I. R. Fisher, 



Z.-X. Shen, PNAS 108, 6878 (2011). 
[4] M. Yoshizawa, R. Kamiya, R. Onodera, Y. Nakanishi, 

K. Kihou, H. Eisaki, and C. H. Lee, Phys. Soc. Jpn. 81, 

024604 (2012). 
[5] T. Goto, R. Kurihara, K. Araki, K. Mitsumoto, M. 

Akatsu, Y. Nemoto, S. Tatematsu, and M. Sato, J. Phys. 

Soc. Jpn. 80, 073702 (2011). 
[6] K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. 

Kontani, and H. Aoki, Phys. Rev. Lett. 101, 087004 

(2008). 



[7] L L Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, [18 

Phys. Rev. Lett. 101, 057003 (2008). 
[8] S. Graser, G. R. Boyd, C. Cao, H.-P. Cheng, P. J. [19 

Hirschfeld, and D. J. Scalapino, Phys. Rev. B 77, 

180514(R) (2008). 
[9] A. V. Chubukov, D. V. Efremov, and I. Eremin, Phys. [20 

Rev. B 78, 134512 (2008). 
[10] H. Kontani and S. Onari, Phys. Rev. Lett. 104, 157001 [21 

(2010). 
[11] H. Kontani, Y. Inoue, T. Saito, Y. Yamakawa and S. [22 

Onari, Sohd State Communications, 152 (2012) 718. 
[12] S. Onari and H. Kontani, Phys. Rev. Lett. 109, 137001 [23 

(2012). 
[13] S. limura, S. Matuishi, H. Sato, T. Hanna, Y. Muraba, [24 

S. W. Kim, J. E. Kim, M. Takata and H. Hosono, Nat. 

Commun., 3, 943 (2012). [25 

[14] N. Fuji war a et aL, unpubhshed. 
[15] T. Nakano, N. Fujiwara, K. Tatsumi, H. Okada, H. Taka- [26 

hashi, Y. Kamihara, M. Hirano and H. Hosono, Phys. 

Rev. B 81, 100510(R) (2010). [27 

[16] M. Sato, Y. Kobayashi, S. C. Lee, H. Takahashi, E. 

Satomi, and Y. Miura, J. Phys. Soc. Jpn. 79, 014710 [28 

(2009); S. C. Lee, E. Satomi, Y. Kobayashi, and M. Sato, 

J. Phys. Soc. Jpn. 79, 023702 (2010). [29 

[17] J. Li, Y. Guo, S. Zhang, S. Yu, Y. Tsujimoto, H. Kontani, 

K. Yamaura, and E. Takayama-Muromachi, Phys. Rev. [30 

B 84, 020513(R) (2011). 



Y. Nakajima, T. Taen, Y. Tsuchiya, T. Tamegai, H. Kita- 

mura, and T. Murakami, Phys. Rev. B 82, 220504 (2010). 

J. Li, Y. Guo, S. Zhang, S. Yu, Y. Tsujimoto, H. Kontani, 

K. Yamaura, and E. Takayama-Muromachi: Phys. Rev. 

B 84, 020513(R) (2011). 

S. Onari and H. Kontani, Phys. Rev. Lett. 103, 177001 

(2009). 

Y. Yamakawa, S. Onari, and H. Kontani, 

arXiv: 1303.2978, 

J. Kunes, R. Arita, P. Wissgott, A. Toschu, H. Ikeda, 

and K. Held, Comput. Phys. Commun. 181, 1888 (2010) 

T. Miyake, K. Nakamura, R. Arita and M. Imada, J. 

Phys. Soc. Jpn. 79, 044705 (2010). 

F. Wang, F. Yang, M. Gao, Z.Y. Lu, T. Xiang, and D.H. 

Lee, Europhy. Lett. 93, 57003 (2011). 

T.A. Maier, S. Graser, P.J. Hirschfeld, and D.J. 

Scalapino, Phys. Rev. B 83, 100515(R) (2011). 

T. Saito, S. Onari, and H. Kontani, Phys. Rev. B 83, 

140512(R) (2011) 

H. Kontani, T. Saito and S. Onari, Phys. Rev. B 84, 

024528 (2011) 

Y. Ohno, M. Tsuchiizu, S. Onari, and Hiroshi Kontani, 

J. Phys. Soc. Jpn. 82, 013707 (2013). 

M. Tsuchiizu, S. Onari, and H. Kontani, 

arXiv: 121 1.3944, 

S. Onari and H. Kontani, unpublished. 



